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+J Abstract 

We have investigated the thermal Hall effect of magnons for various ferromagnetic insulators. 
For pyrochlore ferromagnetic insulators LU2V2O7, H02V2O7, and I^M^Cv, finite thermal Hall 
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conductivities have been observed below the Curie temperature Tq. From the temperature and 
magnetic field dependences, it is concluded that magnons are responsible for the thermal Hall effect. 
The Hall effect of magnons can be well explained by the theory based on the Berry curvature in 
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momentum space induced by the Dzyaloshinskii-Moriya (DM) interaction. The analysis has been 
extended to the transition metal (TM) oxides with perovskite structure. The thermal Hall signal 



was absent or far smaller in La2NiMn06 and YTiOs, which have the distorted perovskite structure 

CN 

with four TM ions in the unit cell. On the other hand, a finite thermal Hall response is discernible 

> 

below Tq in another ferromagentic perovskite oxide BiMnC>3, which shows orbital ordering with a 
larger unit cell. The presence or absence of the thermal Hall effect in insulating pyrochlore and 
perovskite systems reflect the geometric and topological aspect of DM-induced magnon Hall effect. 

PACS numbers: 72.20.-i, 75.47.-m, 75.76.+j 
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I. INTRODUCTION 



The Hall effect is the induction of transverse electric current upon the application of the 
longitudinal electric field. While the Hall effect in nonmagnetic metals or semiconductors is 
usually induced by Lorentz force proportional to the magnetic field, in ferromagnets there 
is an additional component induced by the spontaneous magnetization, which is termed 
anomalous Hall effect.^ Recent theoretical and experimental studies show that the theory 
based on the Berry phase in momentum space accounts well for the observed anomalous Hall 
effects. 1 The spin-orbit interaction gives rise to the topological structure of the Bloch wave 
around band crossing points denoted as magnetic (anti-)monopoles, which act as sources or 
sinks of fictitious magnetic field in momentum space. The anomalous velocity caused by 
the fictitious field is the origin of the anomalous Hall effect. Because the Berry phase is not 
restricted to electrons, the Berry phase induced Hall effect is expected for other particles 
even without the charge. In fact, Hall effects of photons and phonons have been reported 
previously'^' In this paper, we study the Hall effect of magnons, which are the quanta of 
magnetic excitation in magnetic materials. 

Because magnons can carry the spin moments less dissipatively than electrons do, the 
magnon spin current seems to be useful for future spintronics. In this context, the new 
functionalities of magnon spin currents have been investigated recentlyP^ The Hall effect 
may also be useful for the control of the magnon spin current. Theories of magnon Hall effect 
were proposed recently' 11 * 12 ' Fujimoto theoretically suggested that the transverse magnon 
spin current can be induced by the application of the longitudinal magnetic field gradient 
in noncoplanar spin structure.^ On the other hand, Katsura et al. showed that the ring 
exchange interaction induces the Hall effect even in the collinear ferromagnet in the case 
of particular lattice such as the Kagome lattice.^ They also show that the magnon Hall 
effect can be observed with use of heat transport measurment, and derive a formula for the 
thermal Hall conductivity due to magnons.^ 

Quite recently, we have succeeded in the experimental observation of the thermal Hall 
response below the Curie temperature Tq in a ferromagnetic insulator with pyrochlore struc- 
ture Lu2V207.^1 We have found that the temperature and magnetic field dependences are 
consistent with the picture of magnon Hall effect. The observed thermal Hall conduc- 
tivity can be explained by the theoretical model based on the Berry phase due to the 
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Dzyaloshinskii-Moriya (DM) interaction. In this paper, we have investigated the ther- 
mal Hall conductivity in various ferromagnetic insulators to further develop the concept of 
magnon Hall effect. We have found that the thermal Hall conductivities are commonly ob- 
served below Tq in ferromagnetic insulators with pyrochlore structure, In2Mn207, H02V2O7 
as well as LU2V2O7. The temperature and magnetic field dependences for I^M^CV and 
H02V2O7 are similar to the previously observed data of LU2V2O7, except for the variation of 
sign, indicating the generality of our observations. On the other hand, we could not observe 
a finite thermal Hall conductivity for the perovskite ferromagnets YTi0 3 and La 2 NiMn0 6 , 
in which the unit cell contains 4 magnetic transition metal (TM) sites, while a finite signal 
is observed for BiMn0 3 with a larger unit cell including 16 Mn sites. We show that these 
observations reflect the geometric and topological aspect of magnon Hall effect caused by 
the Berry phase due to the DM interaction. 

The organization of the rest of the paper is as follows. In Sec. II, we show the details of 
the sample preparation and the transport measurement. In Sections III and IV, we discuss 
the thermal Hall conductivity caused by magnon Hall effect in pyrochlore and perovskite 
ferromagnets, respectively. Finally, we conclude with a summary in Section V. 

II. EXPERIMENT 

Crystals of LU2V2O7 and YTi0 3 were prepared by the floating zone method. The atmo- 
spheres were Ar gas for Lu 2 V 2 7 and the mixture of Ar and H 2 gases with the ratio of 96:4 
for YTi0 3 . The growth rate was 2 mm/h for both the cases. We prepared polycrystalline 
H02V2O7, In 2 Mn 2 07, and BiMn0 3 samples by high-pressure synthesis. The mixed powder 
of starting material (Ho 2 3 and V2O4 for H02V2O7, ln 2 3 and Mn02 for In 2 Mn 2 07, and 
Bi 2 3 and Mn 2 3 for BiMn0 3 ) with the prescribed ratios was packed into platinum or gold 
capsules (~4 mm0x 6 mm), and heated in a cubic anvil- type apparatus for 60 min. Samples 
were synthesized at 1300 °C and 6.5 GPa for Ho 2 V 2 7 , at 850 °C and 3 GPa for In 2 Mn 2 7 , 
and at 700 °C and 6.5 GPa for BiMn0 3 . The capsules were cooled down to room temper- 
ature before releasing the pressure. Polycrystalline La 2 NiMn0 6 samples were prepared by 
solid-state reaction. Powders of La20 3 , Mn 2 3 and NiO were ground altogether with the 
stoichiometric ratio. The mixed powders were then pressed into half-inch-diameter pellets 
with 2-3 mm thickness and sintered in flowing air at 1100 °C for 48 hours. By means of 
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powder X-ray diffraction, we have confirmed that these samples are without any extra phase 
except for the BiMnOs sample, in which unindexed peaks due to some impurity phase were 
observed but the volume fraction was estimated to be less than 3%. The LU2V2O7 sample 
was confirmed to be a single crystal with use of Laue X-ray diffraction. As for YTi03, we 
have obtained the clear Laue pattern corresponding to the pseudo-cubic perovskite struc- 
ture but the orthorhombic a, b, and c-axis could not be distinguished presumably due to the 
heavily twined structures. 

The magnetization was measured in a Magnetic Property Measurement System (Quan- 
tum Design Inc.). The resistivity was measured with use of Physical Property Measurement 
System (Quantum Design Inc.). We employ a conventional steady-state method for the 
measurements of longitudinal and transverse thermal conductivities. The longitudinal and 
transverse temperature gradients ^ and ^ were measured using both the type E thermal 
couple (T > 20 K) and CX-1050 thermometers (T < 50 K). The longitudinal thermal con- 
ductivity k xx and thermal Hall conductivity n xy were obtained by the following relationships; 

^xx ^ Jq /i\ 

^xx + "Zy U xx §^ 

K 2 ^ 

_ U xy _ UJ X y _ ^xxdy /m 

— . .2 , , .2 ~ . .2 ~ „• ' 



"xx ' ™xy "xx J y 

where u xx , u xy and j g are the longitudinal and Hall components of thermal resistivity and 
thermal current density, respectively. While the transverse temperature gradient |p should 
be antisymmetric with respect to the magnetic field, the if-symmetric component is ob- 
served due to the small asymmetry of the thermal probes in the actual measurements as 
in the transverse voltage of usual electrical Hall measurements. Following the convention 
of electrical Hall measurement, we have subtracted the if -symmetric component using the 
equation = AT (+ H )-^ T (~ H ) ; w here AT and d are the observed temperature difference 
and distance of the transverse thermal probes, respectively. While the Hall conductivity 
at negative fields obtained by this method is merely the copy of the positive field data, we 
plot the data at negative fields just for the clarity of figure. As shown in Eq. (|2]), the small 
value of k xx is desirable in order to estimate the n xy value precisely from the measurement 
of the transverse temperature gradient The present samples certainly have such small 
K rr values as discussed later. 
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III. MAGNON HALL EFFECT IN PYROCHLORE FERROMAGNETS 



In this section, we discuss the thermal Hall conductivity in pyrochlore ferromagnets 
LU2V2O7, H02V2O7, and In 2 Mn 2 07. The results for LU2V2O7 have been published in a 
short paperP^ Here, the detailed data and analyses for all the three pyrochlore ferromagnets 
are presented and discussed comprehensively. Figure [ija) shows the crystal structure of 
pyrochlore oxide A2B2O7. In this figure, the oxygen ions are omitted for simplicity. The 
A and B sublattices are identical with each other. The B sublattice is displaced by half 
a unit cell from the A one. The sublattice structure, commonly called pyrochlore lattice, 
is composed of corner-sharing tetrahedra, and can be viewed as a stacking of alternating 
Kagome and triangular lattices along the [111] direction. 

LU2V2O7 and H02V2O7 are ferromagnetic Mott insulators with one 3d electron per vana- 
dium site. For these materials, the resistivities increase rapidly with decreasing temperature 
and the spontaneous magnetization emerges below the Curie temperature Tq « 70 K as 
shown in Figs. [2ja) and (b), respectively. Spin-polarized neutron diffraction suggests that 
the 3d orbitals are ordered so that they all point to the center of mass of V tetrahedron.^ 
In the orbital ordered state, the virtual hopping process to the higher-lying states stabilizes 
the ferromagnetic order of S = 1/2 magnetic moments.^ For LU2V2O7, the magnetization 
saturates at low field and the saturated value per V atom nearly coincides with 1 Bohr 
magneton (/is), being consistent with the d 1 electronic configuration as shown in Fig. [2^c). 
While only the vanadium ions are magnetic in LU2V2O7, the Ho ions are also magnetic and 
behave as Ising spins in H02V2O7. Therefore, the saturation field is relatively high and the 
saturated magnetization with the component of Ho moment is much higher than 1 /zb in 
H02V2O7. 

In 2 Mn 2 07 is also a Mott insulator J-^ The resisivity is too high to be measured even at 
300 K. There are three 3d electrons per manganese site. The localized S = 3/2 magnetic 
moments are ferromagnetically coupled to each other due to the strong hybridization among 
the In 5s, O 2p, and Mn 3d orbitals.^ As shown in Figs. [2] (b) and (c), the spontaneous 
magnetization is observed below Tc=130 K and the saturated magnetic moment nearly 
coincides with the expected value of 3 /zb- 

We show the longitudinal thermal conductivities for LU2V2O7, H02V2O7, and In 2 Mn 2 07 
in Fig. [27d). According to the Wiedemann- Franz law, the electric contribution of thermal 
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conductivity is less than 10 ~ 4 W / Km below 150 K even for the least resistive sample of 
LU2V2O7, indicating that the heat current is not carried by electronic carriers but by charge 
neutral excitations such as phonons and magnons for these samples in this temperature 
region. The thermal conductivities for these samples are quite small, 2 — 3 W/Km even at 
300 K and decreases with decreasing temperature. These behaviors are frequently observed 
in transition-metal oxides, in which both spin and orbital degrees of freedom are active 
because of the strong electron correlation.^ These are likely because the mean free paths of 
phonons and magnons are suppressed due to the fluctuations of spin and orbital in terms of 
exchange-striction and/or Jahn- Teller coupling. As mentioned above, the small longitudinal 
thermal conductivity is favorable to the sensitive measurement of n xy . 

We reproduce the thermal Hall conductivity for LU2V2O7 in Fig. pp2 A finite thermal 
Hall conductivity is observed below Tq=70 K while it is negligible at 80 K. The magnitude 
shows a maximum around 50 K and decreases with decreasing temperature from 50 K. 
As for the magnetic field dependence, the thermal Hall conductivity steeply increases and 
saturates in the low magnetic field region similarly to the magnetization, which indicates 
that the observed thermal Hall effect is not normal Hall effect proportional to the magnetic 
field but the anomalous Hall effect depending on the direction of magnetization. At low 
temperatures, the thermal Hall conductivity gradually decreases with magnetic field after 
the saturation. 

We now discuss the possible carriers of the thermal Hall current, i.e., either phonons 
or magnons. The thermal Hall effect caused by phonons was previously reported in 
Tb3Ga50i2,^ and theories based on spin-phonon coupling were proposed to explain itP^ 
Nevertheless, we have ascribed the thermal Hall conductivity in LU2V2O7 to the magnon 
Hall effect for the following reasons. Since the mean free path of phonons is expected to in- 
crease with magnetic field due to the reduction of spin-phonon scattering, the decrease of the 
thermal Hall conductivity in the high field region at low temperature cannot be explained in 
terms of the phonon mechanism. On the other hand, the reduction of the magnon popula- 
tion caused by the gap opening due to Zeeman effect will diminish the magnon contribution 
of the thermal Hall conductivity in the high magnetic field. In the theory of phonon Hall 
effect based on the spin-phonon scattering,^ the thermal Hall angle k xv /k xx is expected to be 
proportional to the magnetization. Nevertheless, this is not the case for LU2V2O7. Figure [4] 
(a) shows the temperature dependence of thermal Hall conductivity at ^qH=0.3 T and 7 T 
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being compared with the magnetization. To be precise, we plot the averaged thermal Hall 
conductivity between 0.1 T and 0.5 T for the 0.3 T value and between 6.5 T and 7.5 T for 
the 7 T value because the observed data are scattered. The error bars are estimated by the 
standard deviation divided by the square root of the averaged data number. The thermal 
Hall conductivity is quite small even at [iqH=7 T at T=80 K while the magnetization is 
more than half of the maximum value at this temperature and magnetic field. The difference 
of the temperature dependences in the high magnetic field can be explained by the magnon 
picture because the magnon propagation is caused by the ferromagnetic interaction, and is 
not valid in the magnetic field induced spin polarized state above T c . Thus, the temperature 
and magnetic field dependences are well understood with the picture of magnon Hall effect. 

In order to further examine the scenario of magnon Hall effect, we have also investigated 
the thermal Hall conductivity for other pyrochlore ferromagnetic insulators H02V2O7 and 
In 2 Mn 2 07. In Fig. |5j we show the magnetic field dependence of thermal Hall conductivity for 
Ho 2 V 2 7 at various temperatures. The positive thermal Hall conductivity is observed below 
the transition temperature Tc=70 K. The magnetic field and temperature dependences are 
quite similar to the case of Lu 2 V 2 07 while the magnitude is slightly smaller. In Fig. |4](b), 
we show the temperature dependence of thermal Hall conductivity and magnetization at 
fioH=0.5 T and 7 T. (Similarly to the previous case, we plot the averaged thermal Hall 
conductivity between 0.2 T and 1 T for the 0.5 T value and between 6.5 T and 7.5 T for the 
7 T value.) The thermal Hall conductivity is quite small above Tc even under high magnetic 
field, similarly to the case of Lu 2 V 2 7 . 

The difference from Lu 2 V 2 07 is the presence of Ho / magnetic moment. We show the 
magnetization curve at 20K for Ho 2 V 2 7 in Fig. gc) to compare with the thermal Hall 
conductivity. At this temperature, the magnetization shows a kink structure around 0.3 T 
but gradually increases with magnetic field even above the kink field. The kink corresponds 
to the saturation of ferromagnetic vanadium moments and the increase in the high magnetic 
field is reflected by the gradual alignment process of paramagnetic Ho moments. On the other 
hand, the thermal Hall conductivity saturates at low magnetic field and does not increase 
with magnetic field after saturation, which indicates that the thermal Hall effect is caused 
only by the ferromagnetic and Heisenberg-like vanadium moments. This is consistent with 
the scenario of magnon Hall effect because the magnon picture is not valid for paramagnetic 
and Ising-like Ho moments. 
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Figure 6 shows the magnetic field dependence of thermal Hall conductivity for In2Mn207 
at various temperatures. A finite thermal Hall conductivity is observed below Tc=130 
K. While the sign is negative and the magnitude is large (n xy ~ —2 x 10 -3 W/Km at 
100 K ) in this case, the temperature and magnetic field dependences are quite similar to 
the previous cases. The decrease of thermal Hall conductivity after the saturation is also 
observed at low temperature. In Fig. |2](c), we show the temperature dependence of thermal 
Hall conductivity and magnetization at fioH=0.3 T and 7 T. (Similarly to the previous 
cases, we plot the averaged thermal Hall conductivity between 0.2 T and 0.5 T as the 0.3 
T value and between 6.5 T and 7.5 T as the 7 T value.) The thermal Hall conductivity is 
fairly suppressed above Tq even under magnetic field also in this case. 

The thermal Hall signal in pyrochlore ferromagnetic insulators can be explained by the 
theory of magnon Hall effect based on the Dzyaloshinskii-Moriya (DM) interaction.^ Since 
the midpoint between any two apices of a tetrahedron is not a center of inversion symmetry 
in the pyrochlore structure, there is a nonzero DM interaction 

H DM = J2 D ir(S i xS j ), (3) 

(id) 

where, Dij and Si are, respectively, the DM vector between the sites i and j, and the TM 
spin moment at the site i. The sum is taken over all pairs of neighboring sites. As shown in 
Fig. ffla), the DM vectors on a single tetrahedron of the pyrochlore lattice are determined 



by Moriya's rule^ttSLU anc [ distributed as 

D13 = -^(-1,1,0), £^=-^(-1,-1,0), (4) 

£>43 = ^(0, -1,1), D 12 = -^(0, -1, -1), (5) 

D u = 4kl,0,l),£>23 = -^=(1,0,-1), (6) 



where 1,2,3 and 4 denote the sites shown in Fig. [7] (a). The DM vector on each bond is 
perpendicular to the bond and parallel to the face of the surrounding cube. Let us see that 
the DM interaction does not disturb the ferromagnetic ordering. To this end, we replace Si 
with ((S) +5Si), where (S) denotes the ordered moment. For a single tetrahedron, the DM 
interaction gives rise to the first order term in the fluctuation SSi as 

4 
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where the superscript Kl denotes the single tetrahedron shown in Fig. [7] (a). However, 
since X^(^i) ^ij = f° r an y h one can confirm SHdm = 0. This implies that the collinear 
ferromagnetic ground state is stable against the DM interaction. 

We now briefly explain how the DM interaction can act on the magnons like an effective 
gauge field and gives rise to the Berry phase effect. The Bloch state of a single magnon with 
momentum k is defined by 

i*> = -^E e * Hl i < >' ( 8 ) 

i 

where \i) is the magnon state, in which the spin state at the position Ri is S — 1 with 
all the other spins being completely aligned with the H direction. The matrix element 
corresponding to the transfer integral of magnons is expressed as 

{il-JSi-Sj + Dij-iSiXS^j) 

= -\<i\J(e-«S+Sj + e«STS+)\j) = (9) 

where J is the nearest-neighbor ferromagnetic exchange, and S ± is the operator that in- 
creases (decreases) the spin component along the direction n = H/\H\. The relation 
Je 1 ^ = J + iDij ■ n determines the complex transfer integral for magnons. Note that the 
component of the DM vector perpendicular to n does not contribute to the spin-wave Hamil- 
tonian up to quadratic order. The phase factor due to the DM interaction can be viewed 
as a "fictitious magnetic flux". It should be noted that the lattice geometry is important 
to avoid cancellation of the effect of phase factor; the inequivalent loops in the unit cell are 
necessary for a finite thermal Hall conductivity as suggested previously.^ There are such 
inequivalent loops in the unit cell in the pyrochlore lattice. Figure \f^b) shows the (111) 
plane cross-section of pyrochlore lattice, which is the Kagome lattice. The out-of-plane 
component of DM vector perpendicular to the (111) plane is also depicted in this figure. 
The Kagome lattice is certainly composed of inequivalent loops, namely the triangles and 
hexagons. While the total magnetic flux in the unit cell is zero, the Berry curvature, i.e. the 
fictitious magnetic flux in momentum space, becomes nonzero owing to the inequivalence of 
the loops, which may induce the thermal Hall effect of magnons. 

We now turn to the quantitative calculation of the thermal Hall effect and its compari- 
son with the experimental results. A formula for the thermal Hall conductivity of magnon 
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systems was first proposed in Refs. [T2][T3l However, it has recently pointed out by Mat- 
sumoto and Murakami that an additional term corresponding to the rotational motions of 
magnons is missing in this formula.^ In the following, we will first present the derivation of 
the previous formulsP^^ and show how this formula relates the thermal Hall conductivity 
to the Berry curvature in momentum space. Then we will present the correct formula based 
on the new formulation, which is also expressed in terms of the Berry curvature.^ The ratio 
of the DM interaction to the ferromagnetic exchange interaction D/J will also be estimated 
using the new formula. 

Let us first provide a brief synopsis of the derivation of the previous formula. We start 
from the spin Hamiltonian consisting of the ferromagnetic exchange interactions, the DM, 
and Zeeman terms: 

H = H + # DM , (10) 

H = - JJ2 S i ■ S j - 9HbH ■ (11) 

(m> * 

We can derive the spin-wave Hamiltonian appying the Holstein-Primakoff transformation, 
which is a slight modification of the calculation of the matrix element Eq. Q presented 
above. The spin-wave Hamiltonian in momentum space reads 

tfsw = X>t W (*0**> (12) 

k 

where ^>(k) = (£>!,&, &2.fe, &3,fc ; &4,fc) T are the boson operators that annihilate magnons, and 



/ c i<t>i?. cos ( k v~ k z") g»</>13 cos(k x +ky) c i^ 14 cos(k z -k x ) \ 
cos(0i 2 ) cos(^i3) cos(</>i 4 ) 



H{k) = -2JS 







gj023 cos(k x +k z ) gi0 2 4 cos(k x —k v ) 
cos(4>23) COs((/>24) 







V 



-2043 cos (ky-\-k z ) 
COS(043) 





+ const., 



(13) 



/ 



where k a = k ■ a (a = x,y,z). The lower triangle of the matrix is understood to be filled so 
that the matrix is hermitian. The phase factors faj in Eq. (13) are uniquely determined by 



the DM vectors and the direction of the magnetic field n = (n x , n y ,n z ) through the relation 



tan< 



[Dij ■ n)/ J (see Supporting Online Material of Ref. [13] for details). 



One can diagonalize the spin-wave Hamiltonian at each k and obtain the eigenvector 
\ijj m (k)) with corresponding eigenvalues u m (k), where m (m — 1,2,3,4) is the band index. 
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Using the Kubo formula, one can calculate the thermal Hall conductivity K a @. In the low- 
temperature region, the dominant contribution comes from the small momentum region of 
the lowest (m = 1) magnon band owing to the Bose nature of magnons and the fact that 
the lowest band is well separated from the other bands near k = 0. Retaining only the 
first-order terms in the DM interaction, the analytic expression for the anomalous thermal 
Hall conductivity due to magnons is obtained as 

A 2 f d 3 k /d^(k) 9Vi(*0' 



d 3 k / frMfc) 

2Ii„^ l(fe)Im \l^ 



dk 



a/3 



//|!7 'l.b 



fc„T % 



(14) 



n 3 / 2 ha"~ ' 2JS ' V 2JS' 

where A = 8JS + 2gp B H, p m (k) = [exp(/3w ?n (fc)) — 1] _1 denotes the Bose distribution 
function, and the integral is over the Brillouin zone (BZ). ipi(k) and u>i(k) are, respectively, 
the eigenvector of %{k) corresponding to the lowest magnon mode and its dispersion. In the 



second line of Eq. (14), lA n (z) = YlkLi i s the polylogarithm, a the lattice constant, and 
- eajSjUjD j (8\/2J) with the totally antisymmetric tensor e a/ 3 7 . We have thus obtained 



the expression for K xy in terms of the Berry curvature in the first line of Eq. (14) as a 



result of the approximation neglecting the upper bands. However, as recently pointed out 
by Matsumoto and Murakami^, the thermal Hall conductivity of magnons is, in general, 
expressed by the Berry curvature without such an approximation. By noting that there is 



an additional correction to Eq. (14) that corresponds to the rotational motions of magons, 



they derived the following formula: 

r d?k 



27 yj 



f fUWl /^n(fc) 

3 c 2 (p n {k))lm^-^— 



dip n (k) 



dkr- 



(15) 



where c 2 (p) = (1+p) flog 
the upper bands, we have 



1+p 
P 



2T 



- (log p) 2 — 2Li 2 (— p). Using the approximation that neglects 



d 3 k 



c 2 (pi(fc))Im 



BZ (2vr)3 " \ dk a 

AklT fknT\ b/2 



dk B 



3n 2 ha \2JS J 
1 

c 2 ' 



t 3 ' 2 dt. 



(16) 



Here we have restored the unit h = k B = a/4 = 1. Note that the third line of Eq. (16) is a 
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dimensionless integral. The temperature and field dependencies of K a p are clearly different 
from those of K a p in Eq. (14). 

In Figs. 8 (a)-(c), we show the fitting of the thermal Hall conductivity at low tempera- 



ture by Eq. (16) for the pyrochlore ferromagnetic insulators. The ferromagnetic exchange 
interaction J is estimated by the specific heat datsP'for LU2V2O7 and the mean field values 
J = k-BTc/4:S(S + 1) are employed for H02V2O7 and In 2 Mn 2 07. The field dependence is 



well reproduced by Eq. (16) as shown in these figures. From the fitting, we estimate the 
ratio of the DM interaction to the ferromagnetic exchange interaction D/J as D/J = —0.38 
for Lu 2 V 2 7 , D/J = -0.07 for Ho 2 V 2 7 and D/J = 0.018 for In 2 Mn 2 7 . It should be 
noted that the estimated D/J values have large error bars because it is quite sensitive to 
the estimate of J. (We have the estimate for \D/J\ ~ 0.007 — 0.035 in the case of h^Mj^Oy 
if the error in the estimate of J is within 30 percent.) Nevertheless, the estimated order of 
magnitude \D/J\ ~ 10 _1 — 10 -2 is still meaningful, which is comparable to those reported 
in TM oxidesP^l The difference of the sign of D/J between In 2 Mn 2 07 and the pyrochlore 
vanadates may be ascribed to the different electronic configuration (d l for LU2V2O7 and 
H02V2O7, and d 3 for In^Mn^Oy) while the accurate estimation of many virtual hopping pro- 
cesses is required to determine D and J. The first-principle band structure calculation is the 
most reliable approach to theoretically estimate the D/J as employed by Xiang et al. for 
Y2V2O7 using DFT+U calculation.^ The detailed comparison of the Dzyaloshinskii-Moriya 
interaction between pyrochlore vanadates and In 2 Mn 2 07 with band calculation is a future 
problem. As exemplified in Fig. 4(c), theoretically calculated thermal Hall conductivity k c x ^ c 
steeply increases with temperature. This is qualitatively consistent with the experimental 
observation in the low temperature region. The thermal Hall signal is within the noise level 
at 10 K but clearly observed around 20 K-30 K for all the pyrochlore ferromagnets as shown 
in Figs. 3, 5, and 6. These may be caused by the steep temperature dependence of magnon 
Hall effect. The weaker increase of experimentally observed thermal Hall conductivity in the 
higher temperature region seems to be caused by the effect of magnon-magnon interaction, 
which is important in this temperature range but neglected in the theory. 

As we have seen, the thermal Hall signal due to the magnon Hall effect is commonly 
observed in the pyrochlore ferromagnetic insulators. The order of magnitude and the tem- 
perature/field dependences are quite similar among the materials. The phenomenon of 
magnon Hall effect seems generic in pyrochlore ferromagnets. 



12 



IV. MAGNON HALL EFFECT IN PEROVSKITE LATTICE SYSTEMS 



In this section, we investigate the magnon Hall effect in ferromagnetic insulators with 
perovskite like crystal structure and clarify the effect of lattice geometry on the magnon 
Hall effect. The crystal structure of perovskite oxide ABO3 with GdFeC>3-type orthorhombic 
distortion is shown in Fig. [I] (b). The crystal structure is composed of corner sharing 
B06 octahedra and interstitial A ions. The BC>6 octahedra are tilted alternatively in the 
orthorhombically distorted GdFeCVtype structure. The rare earth and transition metal 
usually occupy the A and B sites, respectively. 

The antiferromagnetic interaction usually works between the magnetic moments of near- 
est neighboring transition metals in perovskite Mott insulators. In order to stabilize the 
ferromagnetic order in the perovskite lattice, the staggered ordering of two different transi- 
tion metals is effective. This is certainly realized in double perovskite oxide L^NiMnOg , in 
which Ni 2+ (t2 9 ep and Mn 4+ (t2 5 ) ions show the staggered alignment in the B-site sublattice.^ 
The crystal structure is monoclinically distorted with the space group of P2\/n in the low 
temperature region. According to the Kanamori-Goodenough rules, the ferromagnetic su- 
perexchange interaction works between Ni and Mn moments. The ferromagnetic sponta- 
neous magnetization is observed below 280 K in L^NiMnOe as shown in Fig. [9]^a). The 
magnetization curve saturates at low magnetic field showing a typical ferromagnetic behav- 
ior (Fig. |9](b)). The saturated magnetic moment is slightly smaller than the expected value 
(2.5 /^B/(Ni,Mn)). The resistivity increases rapidly with decreasing temperature (Fig. [9^c)). 

Another way to stabilize the ferromagnetic order in perovskite oxides is orbital ordering. 
The orbital-order induced ferromagnetic state is realized in YTiOs and BiMnOs. YTiOs has 
one 3d electron per Ti site, which occupies one of the triply degenerate t2 g states. The degen- 
eracy is lifted by the four sublattice orbital ordering with d-type Jahn- Teller distortion) 27 * 28 ! 
We show the temperature dependence of magnetization and the magnetization curve at 5 K 
for YTi03 in Figs. [9] (a) and (b), respectively. The ferromagnetic order is stabilized by the 
orbital ordering below Tc=30 K. The saturated magnetization almost coincides with the ex- 
pected value for the d l electronic configuration (1 /ie/Ti). A neutron diffraction experiment 
observed the gapless magnon spectrum in this material.^ The resistivity is too large to be 
measured even at 300 K, thus not shown in Fig 9(c). 

In BiMn03, four 3c? electrons at the Mn site shows the high spin configuration due to 
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the large Hund coupling. While three t2 9 states with the same spin direction are fully 
occupied, there is only one electron in the doubly degenerated e g states. As a result, the e g 
orbital is ordered with 16 Mn sites in the unit cell as shown in Fig. [lO~P The temperature 
dependence of magnetization and the magnetization curve are shown in Figs. [9] (a) and (b), 
respectively. The ferromagnetic order is stabilized below around 100 K. 131 1 32 1 The saturated 
magnetic moment roughly coincides with the expected value (4 /ie/Mn). The resistivity is 
too high to be measured also in this case. 

Figure [9] (d) shows the longitudinal thermal conductivity for L^NiMnOg, YTiC>3, and 
BiMnCV The magnitude of thermal conductivity is small also for these perovskite ferro- 
magnetic insulators. For L^NiMnOe and YTiC>3, the thermal conductivity monotonically 
decreases with decreasing temperature. On the other hand, it has a broad peak structure 
around 50 K for BiMnC>3. 



We show the thermal Hall conductivity for L^NiMnOg and YTiC>3 in Fig. 11, and for 



BiMnC>3 in Fig. 12 Finite thermal Hall conductivity is not discernible for L^NiMnOg and 
YTiC>3 even below Tq. On the other hand, we have observed the negative thermal Hall signal 
below Tq in BiMn03 . Similarly to the cases of pyrochlore ferromagnets, the thermal Hall 
conductivity for BiMnOs is nearly proportional to the magnetization but tends to decrease 
with magnetic field after the saturation at low temperature. The temperature dependences 
of thermal Hall conductivity and magnetization at HqH=1 T and 7 T for BiMn03 are plotted 
in Fig. [13} (Similarly to the previous cases, we plot the averaged thermal Hall conductivity 
between 0.5 T and 1.5 T as the 1 T value and between 6.5 T and 7.5 T as the 7 T value.) 
The thermal Hall conductivity is suppressed above Tq ~ 100 K even at 7 T while the 
magnetization is large just above Tq in the high magnetic field. Similar to the pyrochlore 
case, these behaviors can be explained in terms of magnon Hall effect. 

The presence or absence of the magnon Hall effect in perovskite oxides can also be well 
explained by the theoretical model based on the DM interaction. In the ideal cubic perovskite 
structure, where the midpoint between two TM ions is a center of inversion symmetry, the 
DM interaction is absent. In many materials, however, the DM interaction is allowed because 
of the distortion, and it is not trivial whether magnon Hall effect can be observed or not. 

As an example, we show the structure and the pattern of DM vectors in GdFeOs type 



distorted perovskite structure in Fig. 14 The a, b, c and x, y, z axes are belong to the 
orthorhombic and pseudocubic coordinate systems, respectively. In this structure, there are 
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four transition metals (1-4 in the figure) in a unit cell and the DM vectors between the 
neighboring sites i and j are distributed as follows: 

D 12 = (a c , (3 C , 0), D u = (-a c , C , 0), (17) 
D 23 = (a ab , -(3 ab , 'j ab ), D 41 = (a ab , -(3 ab , -7 a6 ). (18) 

As described in the previous section, only the component of the DM vectors parallel to the 
magnetic field contributes to the Berry curvature; hence we illustrate the component of DM 



vectors perpendicular to the lattice plane in Fig. 14 (b) (y component in the zx plane). 

In the following discussion, we consider the situation in which the magnetic field is applied 
along the psudocubic y axis and magnon transfers in the zx-plane. When the magnon moves 



on the loop 1— > 2 — > 3 — > A — > I in the Fig, 14 (b), magnon gains the phase factor (0 in the 



counterclockwise direction), whereas a phase factor is — (0 in the clockwise direction) when 
magnon moves on the next neighboring loop 1 — ?• 4 — 3 — 2 — ?- 1 because Dij = —Dji- 
Therefore, the flux, which is the phase factor of magnons, are staggered in a zx plane. In this 
case, the absence of the magnon Hall effect can be roughly understood as follows. Suppose 



that the system is purely two-dimensional as described in Fig. 14 (b). Then if we translate 
the system by a half length of the lattice constant in the x direction and apply the symmetry 
operation that rotates the plane by angle 7r about the x axis, the flux pattern returns to the 
original one. This immediately implies k zx = —k zx and thus the thermal Hall conductivity 
of magnons is zero. Note that the ferromagnetic spin configuration itself is reversed by the 
above it rotation. A similar argument applies to the case where the magnetic field is along 
the x direction or the z direction. A more rigorous justification based on the symmetry of 
the (three-dimensional) spin-wave Hamiltonian is given below. 

In the following, we explicitly show the cancellation by deriving the effective spin-wave 
Hamiltonian for the case of staggered flux pattern and showing that the Berry curvature at 
any fc-point is exactly zero because of symmetry reasons. Let us start from the original spin 
Hamiltonian defined by 

H = ^[JijSi ■ + ■ (5, x S,)] - 9f i B H ^ Sf, (19) 

<y> * 

where denote nearest neighbor pairs. Here, we consider the ferromagnetic exchanges 
with Jij = —J a b on the xy-pl&ne and those with = — J c along the z-axis. Note that 
we have neglected further neighbor exchange interactions. We then take the y-direction 
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as a quantization axis of spins, and apply the standard Holstein-Primakoff transformation, 
yielding the spin-wave Hamiltonian as 



D(k) 
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Dt(fe) 



+ const. 



(20) 
(21) 



where = {bi )k M,kM,kM,k) T and 

D(k): 



J c e 1 ^ cos fc z J a fe(cos k x + cos 



(22) 



J a b(cos k x + cos fcy) J c e** cos /c 2 

with = arctan(/3 c / J c ) — arctan(/3 a {,/J a ;,). Here we have performed a gauge transformation 
which makes the Hamiltonian less cumbersome. Note that k a = k ■ a (a = x,y,z), where a 
corresponds to half the lattice translation in the a-direction. 

We are now ready to show that the Berry curvature for each band at any k is exactly zero. 
To this end, we make use of the symmetry of the spin-wave Hamiltonian. The Hamiltonian 
%(k) is invariant under the following transformation: 



(£ x H(fc)£T = ft(*0, 



(23) 



where the matrix S x is defined by 



(24) 



/o o i o\ 

1 

10 

\0 1 0/ 

and asterisk (*) denotes the complex conjugate. Note that this symmetry is not present 



in the spin- wave Hamiltonian for the pyrochlore ferromagnet Eq. (13). If we suppose that 



the eigenvalue of "H(fc) is nondegenerate, we can determine the form of the eigenvector that 



should be invariant under the transformation in Eq. (23) 



\ip m (k)) = [u m (k),v m (k),u* m (k),v* m (k)] r - 



(25) 



where m (m = 1,2,3,4) are the band indices. For the mth band, the Berry curvature is 
given by 



-21m / %® 
\ ok a 



d^ m (k) 



dk f 



(26) 
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However, due to the special form of the eigenvector Eq. (25), one can easily see that the 
right-hand side of Eq. (26) is always zero. This proves that F^\k) = for any k. The 
same symmetry applies to the case where the magnetic field is applied along the z-axis or 
x-axis. Thus, we conclude that the thermal Hall conductivity n a p should vanish in this 
system irrespective of the direction of the field. 

YTiC>3 certainly has the GdFe03-type orthorhombic crystal structure. While the crystal 
structure of L^NiMnOg is monoclinic, it can be approximately viewed as GdFeOa structure 
because the difference of monoclinic angle from 90 degree is less than 0.1 degree .^Therefore, 
the absence of thermal Hall effect in YTiC>3 and L^NiMnOg is well explained by the above 
theory. On the other hand, BiMnC>3 have a larger unit cell with 16 Mn sites. In this case, 
the cancellation of the Berry curvature due to symmetry reasons may be avoided, and thus 
a finite thermal Hall conductivity is expected and actually observed. 



V. SUMMARY 



In summary, we have investigated the magnon-induced thermal Hall conductivity in var- 
ious ferromagnetic insulators. In addition to the previously reported LU2V2O7 case™ the 
finite thermal Hall conductivity is also observed in other pyrochlore ferromagnets H02V2O7 
and In 2 Mn 2 07. The temperature and magnetic field dependencies as well as the order of 
magnitude are quite similar to those of LU2V2O7, which indicates that the observed thermal 
Hall conductivity in the pyrochlore ferromagnets can be generically ascribed to the magnon 
Hall effect due to Dzyaloshinskii-Moriya interaction. In the perovskite ferromagnets, the 
thermal Hall signal is indiscernible for L^NiMnOg and YTiC>3 , which contains 4 transition 
metal(TM) sites in a unit cell, but a finite signal is observed for BiMn03 with the larger 
unit cell with 16 TM sites. The presence or absence of magnon Hall effect can also be 
well explained by the theory based on the Dzyaloshinskii-Moriya interaction. Our study 
further revealed that the lattice geometry affects the topological strucuture of the Bloch 
wavefunction and is thus essential for the observed magnon Hall effect. 
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Figures 



(a) Pyrochlore oxide A 2 B 2 7 (b) Perovskite oxide AB0 3 




FIG. 1: (Color online) Crystal structures of (a) pyrochlore oxides A2B2O7 and (b) perovskite oxides 
ABO3. The oxygen ions in (a) are omitted for simplicity. 
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FIG. 2: (Color online) (a)-(d) Magnetic, electric, and thermal transport properties for LU2V2O7, 
H02V2O7, and In2Mn207. (a) Temperature dependence of the resistivity, (b) Temperature depen- 
dence of the spontaneous magnetization (fioH = 0.1 T). (c) Magnetization curves at T = 5 K for 
all the samples. For H02V2O7, the 20 K data are also shown, (d) Temperature variation of thermal 
conductivity. 
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FIG. 3: (Color online) Thermal Hall conductivity as a function of the magnetic field for LU2V2O7 
at various temperatures. The magnetic field is applied to [100] direction. 
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FIG. 4: (Color online) Temperature dependence of the thermal Hall conductivity and magnetization 
(a) at /iqH = 0.3 T and [XqH = 7 T for LU2V2O7, (b) at hqH = 0.5 T and [iqH = 7 T for H02V2O7, 
and (c) at ^lqH = 0.3 T and [iqH = 7 T for In.2Mn.2O7- K^j c is the theoretically calculated thermal 
Hall conductivity for I^M^Cv at 0.3 T. 
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FIG. 5: (Color online) Magnetic field variation of the thermal Hall conductivity for H02V2O7 at 
various temperatures. 
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FIG. 6: (Color online) Magnetic field variation of thermal Hall conductivity for In2Mn207 at 
various temperatures. 
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FIG. 7: (Color online) (a) DM vectors in pyrochlore lattice and (b) magnetic flux due to the DM 
interaction in the (111) plane of the pyrochlore lattice (Kagome lattice). 




FIG. 8: (Color online) Magnetic field variation of thermal Hall conductivity (a) at 20 K for LU2V2O7 
(b) at 20 K for H02V2O7 and (c) at 33 K for In2Mn207. The solid line indicates the magnetic field 
dependence given by the theoretical formula based on the Dzyaloshinskii-Moriya interaction (Eq. 



(16)). 





finH= 0.1 T 






La 3 NiMnO s 






BiMn0 3 \ 



-(c) 



100 200 
T(K) 



§10* 

s 10 2 
a. 10 1 - 



1C U , 



100 200 

r(K) 




La 2 NiMn0 6 H 
100 200 300 

r(K) 



23 



FIG. 9: (Color online), (a)-(d) Magnetic, electric, and thermal properties of ferromagnetic per- 
ovskite oxides L^NiMnOg, BiMnOs, and YTiO"3. (a) Temperature dependence of magnetization 
at a magnetic field [1qH = 0.1 T. (b) Magnetization curves at T = 5 K. (c) Temperature variation 
of resistivity, (d) Temperature variation of longitudinal thermal conductivity. 




FIG. 10: (Color online). Orbital ordering in BiMnOs. 
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FIG. 11: (Color online). (a),(b) Magnetic field variation of the thermal Hall conductivity for 
(a)La2NiMn06 and (b)YTi03. For YTiOs, magnetic field is applied to the [100] direction of 
pseudocubic crystal. 
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FIG. 12: (Color online). Magnetic field variation of the thermal Hall conductivity at various 
temperatures for BiMnO"3. 




FIG. 13: (Color online). Temperature dependence of the thermal Hall conductivity and magneti- 
zation at [iqH = 1 T and [iqH = 7 T for BiMnOs. 
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(a) 




FIG. 14: (Color online), (a) DM vectors in perovskite structure with GdFeO"3-type distortion. For 
notational convenience, the eqivalent sites 3 and 5 (4 and 6) are distinguished. The a, b, c and 
x, y, z axes belong to the orthorhombic and pseudocubic coordinate systems, respectively, (b) DM 
vectors and flux pattern in the pseudocubic zx plane. The positive direction of flux is taken to be 
counterclockwise. 
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